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There exists a large number of experimental and theoretical results supporting the picture of 
macroscopic qubits implemented by nanoscopic Josephson junctions of three different types - charge 
qubit, flux qubit and phase qubit. The standard unified description of such systems is based on the 
formal quantization of the phenomenological Kirchhoff equations for the corresponding circuits. In 
this paper a simplified version of the BCS theory for superconductors is used to derive microscopic 
models for all types of small Josephson junctions. For these models the state-dependent individual 
tunneling of Cooper pairs couples ground pair states with excited pair states what leads to a more 
complicated structure of the lowest lying energy levels. In particular, the highly degenerate levels 
emerge, which act as probability sinks for the qubit. These models allow also for the coupling to 
phonons as an efficient mechanism of relaxation for all types of junctions. The alternative formulas 
concerning basic spectral parameters of superconducting qubits are presented and compared with 
the experimental data. Finally, the question whether small Josephson junctions can be treated as 
macroscopic quantum systems is briefly discussed. 
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I. INTRODUCTION 



In the last decade remarkable experiments were performed involving measurements and manipulations of states 
for a single or several nanoscopic Josephson junctions (JJ) which were consistently interpreted in terms of two level 
quantum systems^"— . The main assumption in the theoretical analysis is that such a many-body mesoscopic system 
can be effectively treated as a quantum system of a single degree of freedom typically described by a large spin 
or nonlinear oscillator model. The standard construction of the quantum Hamiltonian for different types of JJs is 
based on the formal quantization of the Kirchhoff equation for the corresponding macroscopic circui1>i. Another 
approach uses an effective picture of Bose-Einstein condensate which represents Cooper pairs at low temperatures^^. 
More fundamental derivations based on path-integral techniques lead to essentially equivalent models as welli^ii^. 
The simple representative of this class of single degree of freedom models is the following Josephson Hamiltonian 
describing the effects of Coulomb repulsion and tunneling for a Cooper pair box (CPB)^ 



Here, J^, fc = 1, 2, 3, are spin operators for the spin-j such that the number of Cooper pairs at equilibrium is close to 
J, Ec is the charging energy, Ej is the Josephson energy characterizing the tunneling magnitude, and Ug is a control 
parameter. 

The obtained nonlinear Hamiltonians yield the structure of two lowest energy levels which at the low enough 
temperatures can be separated from the others to form an effective macroscopic qubit. In particular for the Hamiltonian 
(H)) and under the condition Ec S> Ej these states are approximatively spanned by the eigenstates |mo) and |mo + 1) 
of the operator J3 with mo < rig < mo + 1. 

The main problem with such models is the presence of a typically strong and collective coupling to an environment. 
Namely, it is expected that the observed states should be rather well-localized semiclassical ones (like coherent states 
for the model ([T])), which seem to be the only relatively stable with respect to external noiseii. However, the 
semiclassical states for a model of small JJs are characterized by large charge fluctuations which are not observed 
in the experiments with CPBs . Therefore, either environmental decoherence producing semiclassical states does 
not work for Josephson qubits at the typical time scale of the experiments or the standard single degree of freedom 
model is not correct. The first alternative seems to be unlikely because the semiclassical character of observed states is 
confirmed in the recent experiments on atomic Bose-Einstein condensate (BEC) in a double- well potential--. Although 
this is a physically different system its mathematical description is the same as for the standard model of a CPB 
and given by the Hamiltonian of the form ([1}. Therefore, we follow the second possibility and propose an essentially 
modified theoretical description of small JJs. 

In Sec. II two basic approximations for the Bardeen-Cooper-Schrieffer (BCS) Hamiltonian: Bogoliubov-Valatin 
model Hamiltonian and the collective spin models are discussed. It is argued, that with the specific choice of model 
parameters different from the standard ones but consistent with phenomenology of superconductivity, the collective 
spin model better describes physics of small superconducting devices. In particular, this model predicts the existence 
of specific excitations called excited pairs which can be produced by the individual tunneling or scattering from the 
ground pair states. One should stress, that usually excited pairs are not considered in the literature on JJs, in 
contrast to the single electron/hole quasi-particles. Moreover, in the standard derivation the collective character of 
scattering/tunneling is also assumed. 

In Sec. HI this new picture is applied to CPB, in Sec. IV to flux qubit (FQ ), and in Sec. V to current biased 
junction (CBJ). The main feature is the appearance of the additional energy levels, among them highly degenerate 
ones. This level structure essentially modifies the qubit picture of small JJs and the dynamics of their relaxation. 
Section VI is devoted to the comparison of predictions of the presented unifying microscopic model with experimental 
data for all types of superconducting qubits. The advantage of the presented theory is that the approximative formulas 
for the qubit frequencies contain less free parameters than the standard ones but nevertheless agree with the data. 



The simplest version of the Hamiltonian that incorporates pairing interaction and reproduces not only the basic 
phenomenology of superconductivityi^ but is successfully applicable to small superconducting grains^ is the following 




(1) 



II. BCS MODEL AND ITS APPROXIMATIONS 



H — Ho + -ffrcd, 



(2) 
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where 



and 



Ho= ~ '^P^ ^k^Cka (3) 



k,a=± 



H,^d = -'^Ycl^cl_ci^ci+. (4) 

k,l 

Here we denote by |fc,±) a single electron basis of pairs of time-reversed states enumerated by a discrete index 
k = 1,2, . . . , K and by Ck±, c^_|_ the associated set of fermionic annihilation and creation operators. Those single 
electron states correspond to eigenstates with the single electron energy levels within a cut-off hujcut around the 
Fermi energy ep. The coupling constant g > has here a dimension of energy^. 

Although there exists a remarkable exact solution for this Hamiltoniar^ it is too complicated for our purposes. In 
the following we briefly discuss and compare two well known approximative schemes. 

A. Bogoliubov-Valatin model Hamiltonian 

The basic idea of this approach is to replace the Hamiltonian ([2|) containing four-body interaction (jl]) by a two-body 
(mean-field) one 



k,a=± 

-Y,(AcU,_+A*ck-c,+) (5) 



k,a=± 
k 

with the gap parameter A determined using a self-consistent averaging 

^ = ;|E<^^-^fc+>-- (6) 

k 

In particular one can put the model Hamiltonian ^ into the expression for the thermal average to obtain the value 
of A as the solution of the gap equation. It can be done because the Hamiltonian ([5]) can be transformed to the form 
of the Bogoliubov-Valatin free Hamiltonian for the quasi-particles 

Hbv = ^ EkCl^Cka, (7) 

fe,cr=± 

where — \J [tk — £f)^ + A^, described by a set of new fermionic annihilation and creation operators C-k± , Cl,j^ 
being linear combinations of the old ones. The gap parameter A can be chosen positive and the gap equation for the 
thermal average in the grand canonical ensemble determines the temperature dependence A(r). One should notice 
that the gap equation makes sense only for the reference states with indefinite number of Cooper pairs, like grand 
canonical ensemble or coherent-type states which display essentially normal fluctuations. 

The fundamental results of this theory are the relations between the critical temperature Tc (A(rc) = 0) and the 
parameters of the model. In terms of the parametrization used here they read^^ 

A(0) = 1.76fcsTe, (8) 
2 sinh hwcut/g' 

Only the first relation has been tested in many experiments and found to be reasonable while the second one strongly 
depends on the unknown parameter LUcut- 

As shown already by Thirring and Wehr l^'^i^^ the predictions based on the model Hamiltonian ([7]) become exact 
only in the thermodynamic limit for the grand canonical ensemble when the fluctuation terms of the order \/K become 
negligible in comparison to the bulk terms proportional to K . Therefore, one cannot expect that this model gives 
precise structure of few lowest lying energy levels in the case of mesoscopic samples with essentially fixed number of 
Cooper pairs. 
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B. The choice of parameters 

In the modern literature one usually identifies the cut-off frequency cocut with the Debye frequency uiet^ {hujD ^ 
ksTc) what due to © implies 

hhJD = fii^cut > .9 > ksTc. (10) 

In the original BCS paper— the choice fnocut — ^bTc was suggested what implied that all relevant energy parameters 
of the model could be of the same order of magnitude 

A(0) ~ keT, ~ huj^ut ^ 9- (11) 

The choice of ujcut determines important, although not directly measurable, quantities. The first one is the actual 
magnitude of the coupling constant g. The second one is the number of Cooper pairs in the sample at zero temperature, 
given by the formula K/2 = hujcutN{Q) where A''(0) is the density of electronic states at the Fermi surface (spin is not 
counted). Hence the first choice ^ yields this number larger by two orders of magnitude than the second one (fTU]) . 
The are few heuristic arguments in favor of the choice (fTTj) : 

1 . It implies that there is a single energy scale determining the superconducting phenomena given by the magnitude 
of the electron-phonon coupling, while the choice ([TOl) involves two quite different energy scales {hioo ^ kBTc). 
In a general case of particle interaction mediated by bosons the cut-off in boson momentum p^ut influences the 
magnitude of interaction for interparticle distances d < h/pcut- In the case of phonons this yields the distances 
d < a {a is the lattice constant). On the other hand, the "size" of a Cooper pair is at least two orders of 
magnitude larger than a and hence the Debye cut-off should not enter the effective interaction between Cooper 
pairs in the BCS Hamiltonian Q. 

2. For a normal state of a metal the approximative number of thermally excited electrons is ~ kBTN{0) while the 
others are "frozen in a Dirac sea". When the temperature decreases to the critical value Tj,, it is plausible to 
expect that only the ^ kBTcN{0) excited electrons begin to form Cooper pairs and the others remain "frozen". 
Finally, when the temperature approaches zero all of them recombine into Cooper pairs. Hence, the number of 
Cooper pairs at zero temperature given by K/2 should be rather of the order of kBTcNifi) than fiLJBN{Q). 



C. Collective spin model of superconductor 

In order to produce a simple, exactly solvable model Anderson^^ and independently Wada et al^ considered a 
simplification of the Hamiltonian ([2]) neglecting the kinetic energy term Hq- 

Introducing the pair operators bk — Ck-Ck+ and = c^+c^_ one can treat the system as hard-core bosons or 
equivalently as a system of K spins-1/2 with spin operators sj?, a = a;, y, z, such that 

bk^sl~isl = s-, 

bl = Sl+^Sl = 4^ (12) 
1 

2' 

Defining the collective spin operators J = (J^, Jy, Jz) and Ja = '^'^'^ '^^'^ write the Hamiltonians 

Ho = ^ 2(efc - €F)h\bk = X! ^(^'^ ^ + const, 

k k 
i?rcd = -^fj'-if + i.) (13) 



blk = si + 



K 

which are equivalent to the Hamiltonians ([2]|4]) only on the states not including excitations in the form of unpaired 
electrons. In the following we assume that — a-nd use Hj-^d (I13p . called strong coupling limit Hamiltonian, as 
the approximative Hamiltonian of the system. 

The -ftT-spins Hilbert space can be decomposed into subspaces corresponding to irreducible representations of SU (2) 
of the dimension 2j + 1 and multiplicity rj represented by the suitable Young frames^S 

K/2 

C^'' = C^J+i^C^. (14) 

3=0(1/2) 
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One can use as eigenvectors of Hied the orthonormal basis \ j, m; r) 

J2|j,m;r) = j(j + l)|j, m; r), 

Jz\j,m;r) ^ m\j,m;r), (15) 
where r = 1, 2, . . . , r^, to obtain the corresponding eigenvalues of -ffrcd 

Ej.m^'j^ijij + l)-mim-l)). (16) 

For a fixed total number of Cooper pairs N = K/2 + m and hence fixed m the single ground state is given by 
a nondegenerate eigenvector \K/2,m). The highly degenerate states with j — K/2 — p, p — 1,2,..., describe p 
excitations caUed excited pairs in the original BCS paper—. They are still composed of Cooper pairs but their wave 
functions possess different symmetry with respect to permutations of pairs (ground state is completely symmetric) 
given by the corresponding Young tables. This simplified model with an additional structure which takes into account 
single electron excitations has been studied at finite temperatures^^. 

D. Validity of the collective spin model 

In the following we argue that for the mesoscopic samples the structure of the lowest lying levels is reasonably 
well-described by the Hamiltonian ([T^ with the kinetic part Hq treated as a "small" perturbation of /frcd- This is 
consistent with the choice (fTTj) and implies the relatively low number of Cooper pairs K/2 ~ N{0)kBTc. 

In the case of a small electrode, when the Coulomb repulsion should be included, the dominating part of the 
Hamiltonian reads 

Hf,i - - - Jz + Jz) + ^Ec (Jz-fh)\ (17) 

Here Ec = 2e^/C is the charging energy related to the capacitance C, and the parameter m, |m| ^ K, determines 
the average excess number of Cooper pairs in the system. In the case of relatively large system, i.e., for Eq ^ 5, 
the lowest lying states are of the form 1^^/2,771) with |m — m| ^ K. Therefore, the relevant Hilbert space can be 
represented by a highest spin Hilbert space of the dimension K -\- 1 and the Hamiltonian has the same form as (jl7p 
with the collective operators restricted to this subspace. The highest spin Hilbert space is invariant under the action 
of collective operators {JaiCt = x,y^z}. The same holds if the external action on the system is described by the 
Hamiltonian being a function of {Jq} what leads to a large spin model of a mesoscopic J J. This model is essentially 
equivalent to all single-degree of freedom models of CPB used in the literature^^. However, for Ec > g or/and 
non-collective interactions with an environment the states with j < K/2 containing excited pairs become important. 

For the illustration consider the ground state and the lowest excited levels of the Hamiltonian (|17l) with Ec ^ ff, 
the even value of K and to = 0. The ground state has form \K/2, 0) with the energy —g{K/A + 1/2) and the first 
excited state is {K — l)-fold degenerate \K/2— 1,0; r), r = 2,...,K, with the energy —g{K/4:— 1/2) and hence 
separated from the ground one by the energy gap g. The Young tables corresponding to the ground state (upper one) 
and to the first excited states labeled by 6 = 2, 3,..., if (lower one) are as follows 

[l I ■ I ■ |a| b| c| ■ I ■ 
1 | ■ I ■ |q|c| ■ I ■ |jf| 

(18) 

To describe the excitations in the form of unpaired electrons one has to build a more complicated Hilbert space being 

a direct sum of the Hilbert spaces C^^ . Here K' = K — 2p corresponds to the different sets of electronic states ±) 
with p pairs excluded by the Pauli blocking effect^. This extended model has been analyzed at finite temperaturesi^ 
and the main predictions can be summarized as 

A(0) = g = 2kBT,. (19) 

In the following we assume also that 

fuvcut = 5 = 2kBTc (20) 



what is consistent with the BCS prediction © 
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The numerical values of the important measurable relations obtained from the collective spin model differ slightly 
from the BCS values. These deviations are due to neglecting the kinetic energy term which, even in the regime ([TT|). is 
still comparable to the pairing energy. Nevertheless, the collective spin model possesses certain features which seem to 
reproduce better the physics of small superconductors at the temperatures much below the critical one and coupled to 
external fields. First of all the eigenstates of the Hamiltonian p7|) are also eigenvectors of the electric charge operator 
which is a well-controlled observable for small electrodes. Moreover, in contrast to the Bogoliubov-Valatin model 
Hamiltonian, here the states containing excited pairs are manifestly separated from the states containing unpaired 
electrons, not by the energy difference, but by a certain selection rule. Indeed, the Cooper pair states are invariant 
with respect to time reversal operation while the states containing unpaired electrons are not. Similarly, the collective 
spin Hilbert space is invariant with respect to all Hamiltonian perturbations depending on Cooper pair ("spin") 
operators sj? only. Among them there are: 

1. Kinetic energy term 

2. "Local" electric potentials of the form 

U = Y,Uksl (21) 

3. Cooper pair tunneling Hamiltonian from/to an external reservoir given by 

1 ^ 

Trcd^^Y.'^^kst +I3ks]:), (3k eC, (22) 
fe=i 

4. Scattering of Cooper pairs or tunneling through an internal junction, governed by 

K 

H.= J2TmHS7- (23) 

Therefore, at low temperature regime when the number of thermally excited unpaired electrons is negligible, and 
under the assumption the physically relevant states of a small superconductor can be well approximated by the 
low energy sector of the considered collective spin model. 



III. COOPER PAIR BOX 

A Cooper pair box called "charge qubit" is a circuit consisting of a small superconducting island with a small 
capacitance C connected via Josephson junction to a large superconducting reservoir. The Hamiltonian of the isolated 
small electrode has form P7)) with the charging energy Ec comparable to kTc and hence by ([TO)) to g. Notice that 
for a typical CPB our assumption (fTTj) gives K ~ 10'* while the standard one ^TU\\ yields K ~ 10®. 

In the following we restrict ourselves to the energy levels corresponding to the two highest eigenvalues of given 
by j(j + 1) with j = K/2 and j = K/2 — 1. The corresponding eigenvectors and eigenvalues have the following form: 
the nondegenerate level \K/2, m) with the energy 

E{K/2, m) = _2iE±2l + JLm{m - 1) + AEcim - mf , (24) 
4 K 

and the [K — l)-degenerate level |if/2 — 1, m; r), r — 2, . . . ,K , with the energy 

E{Kl2~\,m)^E{Kl2,Tn)+g. (25) 

Figure 2 presents the energies E{K/2, m) and E{K/2 — 1, m) for m = —1, 0, 1 and K » 1 as functions of the control 
parameter m in terms of the energy scale Ec and the dimensionless parameter g = gj^Ec = 0.4 which is close to 
typical values for most of the CPB's implementations. In this limit all levels are equidistant. 

Due to the Coulomb blockade described by the second term in (jl7p one can consider only the states with m = mo, 
Too + 1 satisfying ffi e [TOo,mo + 1]. Taking into account that toq « K and extracting the irrelevant common 
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m 



FIG. 1. The energies E{K/2, m) + g{K + 2)/4 (solid line) and E(K/2 - 1, m) + g{K + 2)/4 (dashed hne) for different values 
oi m, K » 1 and g — 0.4 as functions of the control parameter rh. 



constant 2Ec + 4i?c'[("^o — ™)'^ + (wq — m)] — gK/A we have an effective Hilbert space spanned by the vectors with 
corresponding energies denoted by simphfied symbols 



where s^r = 2, . . . , K and 



|0) = |ir/2,mo), Eo = -2Ec{l-2ng), 

|1) = |X/2,mo + 1), ^1 = 2Ec[l - 2ng - 4.g], 

|s; 0) = \K/2 - 1, mo; s), Wo = -2Ec[l - 2ng - 4g], 

|r; 1) = \K/2 - 1, mo + 1; r), l^i = 2Ec{l - 2ng), 



Ug — m — mo — g. 



(26) 



(27) 



The effective Hamiltonian of CPB 



One should now include the coherent tunneling process through the junction between a small island and a large 
superconducting grounded electrode given by and treated as a small perturbation of PT)) . Introducing the total 
amplitude /3 = X^fcLi Pk one can decompose the tunneling Hamiltonian into collective and individual parts 



J red— -trcd ' -'rod 



1 ^ 

k=l 



X+lm{l3)Jy 



(28) 



The collective part f^^^ preserves the subspaces of the given j. Therefore, if the collective part dominates one could 
consider only the states with j = K/2 to obtain the standard large spin model. To compare the magnitude of the 
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collective component of the Josephson energy E'j = |/3| with its individual counterpart given by 



■ K 

E 

.k=l 



1/2 



K 



n 1/2 



Lfc=l 



(29) 



one can consider a simple toy model with j3k = Ae 



i\k 



0,1,, 



E'j - 1^1 



II 



iXKl 



1-e 

\a\Vk 



< 



,K-1. Then 
2\A\ 



1 



oiA I 



(30) 



what implies for a generic A that Ej ~ Ej\/K 



On the other hand for purely random amplitudes (3k 



|/3fcp and therefore Ej 



Ej. The real system 

should be placed between these two extremal cases of strong interference and random behavior what implies that the 
ratio Ej/Ej increases as a certain positive power of K leading to the domination of the individual coupling. For 
large JJs with small Ec this effect is suppressed by the fact that the level splitting for a fixed j determined by the 
Coulomb repulsion is much smaller than the level splitting for different values of j given by the superconducting gap. 
For small junctions those energy scales are comparable and the individual tunneling prevails. This implies that the 
matrix elements of Tied between the vectors with the same j are negligible in comparison with the elements between 
vectors with \j — j'| = 1. Hence the only relevant matrix elements are the following 



(r;l|T,ed|0) = (r;l|r;,d|0)=??,, 
(s;0|f,ed|l) = (s;0|f;,d|l)=^'. 



(31) 



Since the states \K/2, m) are totally symmetric the matrix elements (l|Sfc |0) (fc = 1, . . . , K) are equal to a constant 
independent of the index k and it is easy to show that (llT^^j^jO) = 0. 

We can write down the full effective Hamiltonian of the CPB including p!7)) and ([51]) which is a direct sum of two 
similar terms 



H, 



CPB 



^CPB 



' ^CPB 



(32) 



acting on the subspaces "H^g and spanned by {|0), |r; 1), r = 2, . . . , K} and {|1), |s; 0), s = 2, . . . , K}, respectively. 
To a large extend both " subsystems" can be treated separately and completely analogically. Therefore, in the following 
we restrict ourselves to the first one with the Hamiltonian (comp. 



E{ng) 



Po. 



f^(|0)(^| 



10 (o|) 



(33) 



Here 



E{ng)^4Ec{l~2ng), Ej 



E 



1/2 



V]t\ 



(34) 



K K 



lO-E^'-i^^i)' ^o = Ei^; (35) 

r=2 r=2 

with = '^'Hr/Ej. The parameter Ej is the Josephson energy describing the transition rate from |0) to the state 
given by a normalized vector |^). The value of Ej can be controlled by external magnetic field and typically Ej < Ec- 
Introducing two vectors |±) 

6 9 

1+) -cos-|e)+sin-|0), 

|-)=cos^|0)-sin^|0, (36) 
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where 9 is defined by cos0 = E(ng) / E{ng)'^ + E^j, and the qubit ohservahles 

'^^ = l+>(+|-|->(-|, (37) 
'^" = l+>(+l + |->(-| 
we obtain from (j33p a new form of the Hamiltonian 

^CPB = \ {^{ng)o^ + E{ng)Po) (38) 

with two eigenvectors |+) and |— ) separated by the energy difference 

w(n,) = ^[iEc{l-2ngW+Ej (39) 

and the third {K — 2)-fold degenerate level corresponding to Pq with the energy E{ng)/2. This third level lies always 
between |+) and |— ) as ~lu/2 < E{ng)/2 < ljj/2. The relevant energy levels as functions of Ug are showed in Fig. 2. 




0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 2. The energies of the states |+) and [ — ) (dotted lines) and the energies given by (|26|l as functions of rig for g — 0.17 and 
S = Ej/AEc = 0.2. 

The external control is performed by changing Ug and applying a microwave radiation. Obviously, if the system 
is completely isolated the qubit space spanned by |+) and |— ) is invariant with respect to the Hamiltonian and 
the external control yielding the usual model of charge qubit. The third, [K — 2)-fold degenerate, energy level 
corresponding to Pq and the lowest energy levels of the second subsystem become important when the coupling to an 
environment is taken into account. 

IV. FLUX QUBIT 

A flux qubit (FQ) is a small superconducting ring interrupted by one or several Josephson junctions. The main 
difference between this system and the CPB is the existence of the additional quantum number /i which accounts for 
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the quantized circular motion of the Cooper pair as a whole. Consider first a collective spin model of such a ring 
without junction treated as a system ol K x {2iijnax + 1) spins-1/2 with spin operators s^^, a — x,y, z, satisfying 



5^ 5'-' 1 



iSklS. 



(40) 



and cyclic permutations of x,y and z. Here, the index k = l,2,...K/2 corresponds to the internal quantum 
numbers characterizing a Cooper pair in its center of motion reference frame. The second quantum number ^ = 
0, ±1, ±2, ±fj.max labels quantized circular states for the Cooper pairs center of motion with respect to the center 
of the ring. The maximal accessible value of is a consequence of the maximal critical current Iq which can flow in 
the ring at zero temperature. This effect can be described by /x-dependence of the coupling constant due to the 
fact that the states with a nonzero angular momenta generate current / proportional to which modifies the gap 
according to the formula 



or 



91^ = 9 



1 



(41) 



(42) 



To justify (|4T|) one should notice that — g/2 is the pairing energy of a single Cooper pair in the center of mass reference 
frame. When the Cooper pair moves with a velocity v, producing a current density I ^ v, its energy increases to 
—g/2 + const x P = —g{I)/2 due to its kinetic energy and the energy of the created magnetic field. When the current 
density reaches the value /q, such that the effective pairing energy g{Io) =0, it is favorable for the system to relax 
from the superconducting state with a current Iq to a normal state with zero current and higher entropy. 
Defining again the collective spin operators = {J^, JJi^ Jp by 



K 



fcM' a = x,y,z, 



and Cooper pair number operators 



K 

E 

k=l 



1 



7^ J 



(43) 



(44) 



one can write the strong coupling limit BCS Hamiltonian for a ring as 



-El ^ A*ext^ 



2 



(45) 



The first term in the Hamiltonian is the standard mean-field BCS pairing Hamiltonian. The second term is 
also of the mean-field type what is the reasonable approximation for the rings with a thickness not larger than the 
penetration depth for a magnetic field (typically ^ 100 nm). In these cases we can ignore spatial variations of the 
current and derive (|45p using the macroscopic relation between the current density j and the vector potential A in a 
superconductor— 



J = 



and the energy of the current in a magnetic field 

£ = - 



TO 



(46) 



(47) 



Here A/" is a density of superconducting electrons, to is an electron mass and A = |A| is assumed constant along the 
loop. Using now the quantization condition for the magnetic fiux that threads the loop of the length € 



A- 



F<S>o, i^ = 0,±l,±2, 



(48) 
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where $o = ^/2e, one obtains the quantized energy 

Replacing F by the flux operator F — X^/^M^m ^-^d adding a shift fie^t = *I'ea;t/*I'o caused by the external magnetic 
flux ^ext one obtains the second term in (j45p with the inductive energy given by 

El = Kh?/'iml^. (50) 

The above choice of flux operator is consistent with the Onsager hypothesis that the flux generated by the circulating 
charge 2e is quantized in the units of <^iy-. 

Notice that we neglect here the Coulomb repulsion term as the charging energy Ec is much smaller than for the 
CPB and the number of Cooper pairs in the system is fixed. 

The physical Hilbert space is spanned by the joint eigenvectors 

J^I-(iM;"M;'"A')-> =iM(> + i)l--(jM;"A';''M)-). 

with the multiplicity and satisfying the condition 

which determines the total number of Cooper pairs in the system. 



A. Ground state and lowest excitations 



Consider first the case fXext = 0. The unique ground state of the Hamiltonian (I45p can be obtained by minimizing 
the energy given by the first term. Indeed, due to the symmetry = g^^ the contribution from the second term 
automatically vanishes for such minimizers. One can easily show that the ground state has a product structure 

|0) = |...(if/2;n^;l)...), (53) 

where ~ tPm ^^^^ the ocupation probabilities minimize the functional X/i.9A'(PM ^ ^Pm)- The numerically 
obtained shape of the probability distribution is presented in Fig. 3. 

The spectral decomposition of the Hamiltonian (|45l) is very rich and contains different types of manifolds of the 
lowest lying excited states. The first type is obtained by the creation of a single excited pair without changing the zero 
magnetic fiux of the ground state, the second one is spanned by the vectors with the fiux quantum number _F = ±1 
and the third corresponds to both types of excitations. The external magnetic flux ^ext = ^^ext^o^ fJ'ext — 1/2, shrinks 
the gap between the ground state and the states with F — 1 what is necessary to separate two qubit states from the 
rest. On the other hand, similarly to the case of CPB individual tunneling processes create excited pairs. Therefore, 
to construct a proper model of FQ we need only the third type of excited states (flux F — 1 and one excited pair) 
represented by the vectors 

\u; r) = I (if/2; 1) . . . {K/2- ~ 1; l){K/2 - 1; + 1; r) . . . {K/2- 1)) (54) 

with V = —Umax: ~ i^-max + Ij ■••7 l^-max and r = 2, 3, K . 



B. The FQ Hamiltonian 



For small rings the condition El 3> g holds and therefore, in order to produce a qubit, we have to switch on the 
external magnetic flux ^ext — '&o/2, (/iext — 1/2). Then the gap between the ground state and the states ([54]) is 
given by 



Eqv = El{1 — 2flext) + gu+i- 



(55) 
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FIG. 3. Ground state occupation probabilities of different circular states. Notice that probabilities vanish far away from the 
borders ensuring stability of the ground state. 



The presence of the junction(s) can be modeled by the generalization of the Hamiltonian 



K 



k,l — l 



(56) 



where T[k, ulljV] is a hermitian matrix of scattering amplitudes. 

Again, similarly to the CPB case, the collective scattering processes which preserve quantum numbers can be 
neglected in comparison with the individual scattering changing j^. Therefore, the relevant matrix elements of the 
perturbation ([55)1 are the following 



{iy;r\Hs\0) = r,, 



(57) 



Similarly, to the CPB case we can write the effective Hamiltonian of the first qubit (the second one is initialized by 
decreasing fi^xt to 1/2) in the lowest order perturbation and projected on the the 2-dimensional qubit space 



Here 



//fq = - (i?(Ae.t) (IO(el - |0)(0|) + Ej (|0)(e| + 10(01)) 



E{fj.ext) = El{1 - 2flext), 



K 



1/2 



(58) 

(59) 
(60) 

(61) 
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K 

10 =I]I]U|i^;0, i.r=2T^ur/Ej. (62) 

In the above we omit the part of the Hamiltonian which describes the "sink states" decoupled from the qubit states. 
Again we observe the structure of levels qualitatively the same as for the CPB (see Fig. 2) with characteristic qubit 
level repulsion caused by a tunneling. 

The final form of the FQ Hamiltonian reads (compare with (|38II39|) ) 

^FQ = ]^u:{jlextW, 
u:{fie.t) = ^[EL{l-2fi,,t)? + Ey (63) 
The qubit is controlled by changing the external flux jlext'^o and applying microwave radiation. 

V. CURRENT-BIASED JUNCTION 

This type of J J, denoted by CBJ, consists of larger superconducting electrodes than in the CPB device what implies 
that the Coulomb energy Ec is much smaller than the gap parameter g. The electrodes are connected to a current 
source which produces a constant but tunable current /. We propose a microscopic model of CBJ using again the 
reduced BCS Hamiltonian (jl3p and emphasizing the role of excited Cooper pairs. In the first approximation we treat 
the system as a closed superconducting circuit with a low value of El which can support a steady current / and the 
junction acting as a perturbation - a scattering center. The starting point is the effective unperturbed Hamiltonian 
similar to (|45p under the assumption that El <C g{I) (see ((4T|l ). The charging energy (neglected in (j45|) ) is also 
small {Ec 9(1)) and the number of Cooper pairs is not fixed what implies that both quantum numbers ^ and 
correspond to certain essentially classical degrees of freedom. The collective coupling of these degrees of freedom to 
an environment produces, by mechanisms mentioned in Sec. I, semiclassical coherent-like states determined by the 
external conditions. Therefore, the ground state of the system can be written as (compare the structure of eigenvectors 

m) 

|0) = ^^0/(M)V'K(n)l-(jM = K/2;n^ = n;r = 1)...), (64) 

^ n 

where the probability amplitudes ^/(/i) and ipxin) display normal fluctuations around mean values ^ I and K/2, 
respectively. In order to construct the lowest excited states we apply a kind of adiabatic approximation fixing the 
semi-classical degrees of freedom and changing the only quantum one related to a number of excited pairs p. Then 
the structure of strongly degenerate excited states is the following 

kp;p) = ^^'PA^A^K{n)\...{j^, = K/2-p;n^ ^ n;rp)...), (65) 

/I n 

where p <^ K and Vp describes degeneracy. 

Remark: The derivation of above is valid for a single JJ in the phase qubit regime. If the device is designed as a 
loop, the loops inductive energy should be taken into account. The additional quantized energy produces different 
initial flux states replacing the ground state (I64p . They yield the different critical currents as the internal loop current 
adds to the external biased one. Therefore, the initialization of the proper flux state must be done before the device 
can be used as a qubilj^. 

Similarly to CPB and FQ we consider a qubit model including only the ground state |0) and the {K — l)-degenerate 
flrst excited states |r; 1) separated by the energy gap g{I). Again the junction acts as a scattering center given by 
the Hamiltonian ((56)) which couples the ground and excited states as in ([57]) . Repeating the analogical construction 
we obtain the following qubit Hamiltonian 

^CBJ - I (oj{I)a' + giI)Po) , (66) 

where 



(67) 
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The two eigenvectors |+) and |— ) of (|66|) are separated by the energy difference a;(/) an the third {K — 2)-fold 
degenerate level with the energy g{I)/2 corresponds to Pq. This third level lies always between |+) and |— ), as 
-uj/2 < g/2 < oj/2 (Fig. 4). 
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l/'o 

FIG. 4. Qubit levels (solid line) and the sink level (dashed line) as functions of biased current I. 

Similarly to the previous cases the external control is performed by tunning the biased current and the coupling to 
microwave radiation. 

VI. COMPARISON WITH STANDARD THEORY AND EXPERIMENTS 

The main predictions of the proposed unified microscopic model of small J Js concern the structure and parametriza- 
tion of their lowest energy levels. The comparison of our model with standard theories and experimental data is 
presented below for all three types of superconducting qubits. There are also another aspects of our approach related 
to dissipation/decoherence processes which will be discussed in the forthcoming paper. 



A. Charge qubits 

The Kirchhoff equation for the CPB is equivalent to the Hamilton equations obtained from the following classical 
Hamiltonian of a fictitious particle^ 

H ^ ~ ^^^^ - Ej cos X. (68) 
2m 

Here the "position" x corresponds to the phase variable, pg ^ is a control parameter and the "mass" is given by 



m=(-) a (69) 
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The harmonic approximation around the potential minimum gives the value of the plasma freqency 

However, for small CPBs one does not use a quantized version of (j68[) but rather a two-level approximation to the 
large spin Hamiltonian ^ 

H = -^[Ec{l- 2ng)a' + Eja^] . (71) 

The Hamiltonian (j7ip is essentially the same as ours p3p and predicts the same qubit frequency as a function of Ec 
and Ej. Therefore, to test our model one should explore the differences. The consequences of the presence of the 
highly degenerate level (probability sink) for the decoherence processes will be discussed in the forthcoming paper. 
Here we refer to Fig. 2 and the formulas (|26p which show that when the control parameter Ug is ramped from zero to 
the value 1/2 the ground state of the initialized qubit is unstable due to the transitions to the level Ei in the range 
(1/2 — g) < rig < 1/2. A similar instability is predicted in the range 1/2 < rig < 1/2 + g for the case of Ug decreasing 
from 1 to 0. 

The instability of the ground state for 0.3 < Ug < 0.7 was reported by Lehnert et alJ^ and attributed to backaction 
generated by currents flowing by the RF-SET device which is used to measure the charge of the CPB. Here we give an 
intrinsic explanation of this effect supported by the data. Indeed, taking experimental^ value AEc/h = 149 GHz and 
puting g/ks = 2Tc = 2.4 K for Al, one obtains g = 0.17 what agrees very well with the observed instability range. 

B. Flux qubits 

The standard description^ is based on the quantization of the following modification of ((68)) 

H=^-Ejcosx + eJ^^^)\ (72) 

where x = 27r$/$o and x^xt — ^n^ext/^o- Here $ is a magnetic flux that threads the loop and ^ext is the external 
flux applied to the loop. In the standard approach one has to compute the spectrum of the quantized version of (|72|) 
using the quartic approximation which leads to a double well potential picture. One should notice that the frequency 
corresponding to a two level approximation depends on three parameters characterizing the device Ec, Ej, El and 
one control parameter ^ext- Applying the prediction of our model (|63p under the assumption \El(1 — 2/ie^t)| >> Ej 
we obtain a very simple formula for the qubit frequency / 

hf ^ \ELil ~ 2t^ext)\ (73) 

depending on the microscopic parameters and the geometry of the sample (see (jSOp ). 

This result can be compared with the experimental date^. The linear dependence given by (I73|) far enough from 
fj-ext = 1/2 is clearly confirmed and the value E^/h ~ 1.5 x 10^ GHz can be extracted from the Fig. 3.(B) presented 
by van der Wal et al}^ The sample is an aluminum 5 /im x 5 /im loop made of 450 nm wide and 80 nm thick lines. 
Using the value (. — 20/i m and the formula ([50]) we obtain the number of superconducting electrons K — 3.3 x 10^. 
This value can be put into the formula 

K = 2nw,utN{Q) = 2hw,utV-^ (3^2^) , (74) 

where V is a volume of the sample and k is a density of electrons (k = 18.06 x lO^^/cm"^ for Al). The substitution 
yields hwcut/2kB = 1-3 K which is close to the critical temperature Tc = 1.2 K. This is a strong support for the basic 
assumption of our model tkUcut — g = Sfc^Tc (see (jlip ) and gives the first estimation of the zero temperature density 
of Cooper pairs K/2V = 2.5 x 10^^/cm'^ in a superconductor obtained directly from the experimental data. 

C. Phase qubits 

For this type of junction one uses a model equivalent to the fictitious particle moving in a washboard potential with 
the Hamiltonian^ being again a modification of (155]) 

H=- Ej (cosx+—x\ (75) 

2m V Iq J 
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The harmonic approximation around the potential minimum gives the value of the plasma frequency 



hC 



1/4 



(76) 



which can be treated as a rough approximation to the qubit frequencyi. More complicated formulas which take into 
account anharmonicity and involve Ej as an additional parameter are also available. On the other hand our formula 
for the phase qubit frequency obtained by combining (|4ip with (|67p and under the assumption g{I) ^ Ej reads 



co{I) c g 



(77) 



It is important that the formula (j77|) involves the microscopic parameter g while the standard expressions depend 
entirely of the macroscopic ones. 

Firstly, one can check roughly the magnitude of the predicted qubit frequencies w(/). Typically, Ej ^ g ~ A9 GHz 
for Al and g ~ 388 GHz for Nb. In all experiments the biased current satisfies 0.85 < I/Iq < 0.99 what leads, using 
dZZl) to a reasonable range of frequencies 0.98 GHz <uj/2tt< 13.6 GHz for Al and 7.7 GHz <uj/2tt< 108 GHz for 
Nb, respectively. 

More detailed comparison is performed for the several examples from the literature. One should notice that 
comparing the experimental data with any theoretical curve describing u}{I) which depends on free parameters is 
difficult and inacurate because the range of the biased current variation is in all experiments a very small fraction of 
the whole interval [0,/o] (see Fig. 5). 




34.00 34.02 34.04 34.06 34.08 34.10 



Ib(^iA) 



FIG. 5. Spectrum of SQUID DSi'^"^ (Nb). (a) The circles represent the experimental dataSA and the line is fitted according to 
((77)1 with Jo = 34.387AtA [g = 2kBTc ~ 388 GHz), b) The fit in the whole range of variation of lo{I)/2tt[GH z]. 

Table I presents the obtained values of /q using the formula (|77)) and the experimental values of uj{I)/2-k for few 
examples of phase qubits. The theoretical values are compared with the reported experimental values of the critical 
current or in the case SQUID DSi^'^ with the value fited to the standard theory. 
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TABLE I. The critical currents for phase qubits realizations 



Junction 


Jo ^iA (exp.) 


Jo /iA (theory) 


SQUID DSi (Nb)^ 


34.275 (fit) 


34.387 (fit) 


current-biased (Al)^ 


13.33 


13.99 


current-biased (Al)^ 


21 


22.4 


SQUID DSsA (Al)^ 


1.26 


1.27 



One should remember that the measurements of the critical current are not very precise because this parameter 
can vary in time for the same sample. The positions of experimental values of w(/) are also quite sensitive to the 
level repulsion phenomena caused by "parasite" two- level systems present in the environment and interacting with a 
qubili^i^. 



VII. CONCLUSIONS 

The presented approach to superconducting qubits differs from the standard one by referring to the microscopic 
Hamiltonian being a simplified version of the BCS Hamiltonian and avoiding a detour via requantization of the 
macroscopic Kirchhoff's equations. As a consequence, in contrast to the standard results, the obtained formulas 
describing the energy spectra depend on the microscopic parameters. Our choice of the microscopic Hamiltonian 
implies also a nonstandard assumption about the cut-off energy scale of the BCS model. We follow the original 
approach where this energy is of the order of superconducting gap, while in the modern literature one chooses the 
Debye energy which is larger by two orders of magnitude. The computation of the number of superconducting electrons 
in the flux qubit device, based on the experimental data, strongly supports this choice. A number of experimental 
results concerning the energy spectra of different types of superconducting qubits is consistent with our model as well. 

Another feature of our model is the importance of individual scattering/tunneling of excited Cooper pairs which 
differs from the standard picture of independent quasiparticles in the Bogoliubov-Valatin scheme. In particular the 
qubit states in our model are spanned by the ground state and the given single excited Cooper pair state. This is 
conceptually very different from the standard picture of macroscopic quantum systems but on the other hand solves 
the puzzle of missing environmental effects which should produce semiclassical behavior. Last but not least, the 
presented model allows new mechanisms of decoherence due to the excited Cooper pair coupling to phonons and 
existence of probability sinks. Those phenomena will be studied in the forthcoming paper. 
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